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Abstract 

An assessment of the turbulence levels present in the Measurement Science and 
Technology (MSTB) branch’s 2 x 6-inch low speed wind tunnel was made using Particle Image 
Velocimetry (PIV), and a turbulence generator consisting of a grid/honeycomb structure. 
Approximately 3000 digital PIV images were captured and analyzed covering an approximate 2 x 
6-inch area along the centerline of the tunnel just beyond the turbulence generator system. 
Custom software for analysis and acquisition was developed for semi-automated digital PIV 
image acquisition and analysis. Comparisons between previously obtained LTA and LV 
turbulence measurements taken in the tunnel are presented. 



MSTB 2 x 6-inch Low Speed Tunnel 
Turbulence Generator Grid/Honeycomb 
PIV Measurements and Analysis 


Introduction 

Particle Image Velocimetry (PIV) measurements were taken as part of the NASA Langley 
Advanced Subsonic Technology (AST) program to better understand the turbulence levels of the 
Measurement Science and Technology (MSTB) low speed wind tunnel for later use in noise 
reduction experiments. A grid/honeycomb structure was used to produce known turbulence 
levels in the tunnel for measurement/probing with the newly developed digital PIV technique 
The goal was to test the feasibility of using digital PIV to make highly accurate turbulence 
measurements, and to assess the turbulence characteristics of the tunnel in general. A comparison 
of previously made LTA and LV turbulence measurements obtained in the tunnel under a similar 
configuration was also made 

Particle Image Velocimetry (PIV) System Description 

Particle Image Velocimetry is an instantaneous 2-D non-intrusive global velocimetry 
measurement technique that has become well established in the fluid mechanics measurement 
community. A basic PIV system involves taking a double exposure recording of typically pm- 
sized seed tracer particles entrained in a flow of interest. 1 Auto- or cross-correlation analysis 
techniques are then used to measure the local displacement levels of the particle image pairs, 
which can be used to gain access to local velocity levels in the flow given the original pulse 
separation timing between the double pulsed light sheets. 

Beyond providing basic information about spatial distributions of velocity in a flow, 
ensemble statistics are often used to measure the spatially resolved mean flow, standard 
deviations, turbulence statistics, and velocity and vorticity correlation levels, among other 
things. 2 ’ 3 In the present application, estimates of 1st and 2nd order turbulence statistics levels 
were made, in addition to ensemble mean and standard deviation levels. 

Tunnel Description and Turbulence Generator System 

The MSTB low speed wind tunnel is a small subsonic open-air induction type (Figure 1). 
Its test section has a cross section of 50x1 50mm, a length of 600mm, and a maximum velocity of 
15m/s. The test section is constructed of high quality glass to provide optical access for the PIV 
light sheets and imaging system. The tunnel was originally designed to provide a stable and highly 
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laminar flow. Seed is introduced at the open inlet structure. The seed in this test was provided by 
a commercial TSI six jet atomizer using mineral oil with a nominal seed size of 1 pm. 

The MSTB wind tunnel had previously been tested for turbulence levels using Laser 
Transit Anemometry (LTA) and Laser Velocimetry (LV) techniques. 4 In those tests, a turbulence 
generator system was used based on a design similar to one used by Xia et. al. 5 in which various 
honeycomb/screen combinations were used to set up and study turbulence. The honeycomb and 
screen combination used in this set of tests consisted of a course screen located 200 mm upstream 
of a honeycomb section of length 127 mm (Figure 2). The course screen was constructed of a 
1mm diameter wire mesh grid at 12.7mm 2 spacing, which gave an open area percentage of 80%. 
The honeycomb consisted of approximately 100 plastic straws cut to length (127mm), and fit 
snugly within the glass test cell through a removable bottom section. This configuration sets up 
an unsteady flow field behind the honeycomb with known turbulence levels. 

PIV Acquisition Description 

The PIV acquisition system used for this set of tests was designed, built, and operated by 
S.M. Bartram and W.M. Humphreys of NASA LaRC (RTG/FMAD/MSTB). It uses two 
frequency-doubled Nd:YAG lasers operating at 532nm, with approximately 600mj energy levels, 
and 10Hz pulse repetition rates. Co-linear light sheets of approximately 1mm thickness by 
100mm width were formed and directed into the test section area (Figure 3). The imaging system 
consisted of a custom polarizing cube image shifter system and a Kodak Megaplus 1.4 digital 
camera system with the camera lens system chosen to provide approximate 1 :2 imaging. The 
entire imaging system was mounted on a precision translation stage to allow imaging different 
regions along the length of the test section without disturbing the light sheet or tunnel orientation 
Digital images were captured and sent to an IBM compatible PC/AT 486/120 computer for 
storage using custom control software which is described in more detail in Appendix A. 

The region and coverage areas studied in this test are depicted in Figure 4. All stations 
were taken along the tunnel centerline both vertically and in the transverse direction of the flow, 
with the initial station covering just past the honeycomb structure to approximately 19mm 
downstream. Each digitally captured image (1320x1035 pixels) covered a 18.86mm x 14.79mm 
area in the test cell. An approximate 30% overlap was provided between stations, which resulted 
in a translation of 13mm between each station. 14 stations were taken in all providing a total 
coverage area of 15.0mm by 182mm. 

PIV Analysis/Validation Processing Description 

A custom and automated PIV analysis system was used in this effort based on a PC/AT 
486/120 platform incorporating an Alacron/i860 array processing system. The system is currently 
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capable of processing 5 autocorrelation displacement measurements per second, and completes an 
entire image field in 1 minute with 50% oversampling. Details of the digital PIV analysis system 
algorithms developed for this effort are provided in Appendix A. A 128x128 pixel correlation 
size was used to provide the needed accuracy levels. This translates to an approximate 1 ,829mm 
spatial resolution capability, and a 0.914mm measurement resolution size. Image sizes ranged 
from 15-50 pm which covered 3-8 pixels on the digital camera. This, along with seed density 
levels of 3-20 image pairs per mm 2 , resulted in correlation signals with SNR of approximately 
50: 1 . A weighted mean centroiding scheme was used, and provided measurement sensitivity 
levels of 0. 1 to 0.8 pixels, which corresponds to a lower limit of approximately 2 to 10 pm. A 
typical autocorrelation processing step is depicted in Figure 5, which shows a raw 128x128 pixel 
digital image area, its autocorrelation pattern, and the resulting scaled vector output. 

An assessment of the image shift calibration stability of the system was made, the results 
of which are depicted in Figure 6. Station-to-station image shift variations were noticed in both 
the axial and transverse component levels, and were most likely due to minor variations in the 
relative mirror alignment of the image shifter due to mechanical vibration or temperature 
fluctuation levels present during the test. An increase in the axial image shift level and decrease in 
the transverse image shift level of +66.89 pm and -98.8 pm were noticed, respectively, from 
station 1 to station 14, with ensemble station standard deviation levels present of 15-29 pm. 

Data validation and post processing efforts also involved using custom and automated 
algorithms based on a PC/AT 486/120 computer platform. Measurement validation rates ranged 
from 10-60% percent based primarily upon the image quality and seed levels available in the raw 
data. A four-tiered validation routine was used to reduce the analyzed data. This involved a 
bandpass validator operation, followed by a local median validator, a 3x3 local mean filter, and 
finally an image shift subtraction operation. The validation process was automated, relying only 
minimally on operator input, and took approximately 1 second per analyzed data file. Following 
data validation, the ensemble mean/standard deviation, and data convergence algorithms were 
applied to the various data sets. 

Results 

More than 3000 digital images were captured and analyzed for this test. This consisted of 
200 data images and 10 image shift reference images at each of the 14 stations. The resulting 
ensemble averaged mean velocity field near the exit of the honeycomb structure is shown in 
Figure 7. Each straw in the honeycomb structure acted like a small jet which resulted in three 
expanding jet structures in the vector field shown in the figure. An entire straw/jet structure can 
be seen along the measurement centerline with 2 additional jet structures partially visible on the 
top and bottom of the field respectively. Data convergence for this initial measurement station is 
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shown in Figure 8 for nine positions in the field, corresponding to the center vector position, and 
8 neighboring vector positions at 5mm vertically and horizontally displaced positions away from 
the center. The ensemble mean is seen to converge in a reasonable fashion with and ensemble 
number above 100-150 data points. The ensemble mean velocity field patterns for all stations are 
depicted in Figures 9 (axial component contour plot underlay) and 10 (transverse component 
contour plot underlay). Ensemble variance levels are depicted in Figures 1 1 and 13, for the axial 
and transverse components, respectively. Three axial cuts were taken along the length of the data 
set as shown in Figures 1 1 and 13, representing data along the center jet centerline and the jet 
edges between the three jets. These results are depicted in Figures 12 and 14, for the axial and 
transverse components, respectively. 

First and second order turbulence statistics calculations were also determined for the 
entire data ensemble. Turbulence intensities for the <u’u’> axial component and <v v > 
transverse component are depicted in Figures 15 and 17. Three axial cuts were taken along the 
length of both of these data sets and are shown in Figures 16 and 18, respectively. Reynold’s 
Stress and out-of-plane mean vorticity calculations were also made, and are depicted in Figures 
19 and 20, respectively. Turbulence intensity and Reynold’s Stress levels are seen to increase 
sharply just beyond the honeycomb structure, peaking at approximately 20mm downstream, and 
then decreasing gradually to a minimum level at 100 mm. Similar trends can be seen in the 
ensemble variance plots in Figures 12 and 14. Vorticity levels and trends show jet-like structures 
out to approximately 25 mm, with vorticity levels minimizing and stabilizing beyond that. 

Standard deviation comparison plots with previously taken Laser Velocimetry (LV) and 
Laser Transit Anemometry (LTA) data 4 are depicted in Figures 21 and 22. The LV and LTA 
data was taken in the same MSTB wind tunnel under similar conditions. Good agreement is seen 
between the different techniques beyond 25 to 30 mm downstream. The PIV data, however, 
showed somewhat higher levels closer to the honeycomb structure. This is most likely due to the 
fact that PIV is a global velocimetry technique and LV and LTA are point measurement 
techniques. As seen in the mean velocity and standard deviation plots of Figures 9, 10, 11, and 
13, the levels are seen to vary considerably in the jet structure regions. The likelihood that the 
PIV axial cut position along the jet edge, and the LTA and LV point measurement locations in the 
transverse directions were coincident, is not very good. General trends and standard deviation 
levels are in good agreement between the three techniques, however. 

Summary 

An assessment of the turbulence levels present in the Measurement Science and 
Technology (MSTB) branch’s 2 x 6-inch low speed wind tunnel was made using Particle Image 
Velocimetry (PIV), and a turbulence generator consisting of a grid/honeycomb structure. 
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Approximately 3000 digital PIV images were captured and analyzed covering an approximate 2 x 
6-inch area along the centerline of the tunnel just beyond the turbulence generator system. 
Custom software for analysis and acquisition was developed for semi-automated digital PIV 
image acquisition and analysis. Comparisons between previously obtained LTA and LV 
turbulence measurements taken in the tunnel are presented. 


Appendix A 

Two software packages were developed for this effort. These included a digital PIV 
(DPIV) acquisition program and digital PIV analysis program for use with a Kodak Megaplus 
vl.4 digital camera system. The software runs under Windows 3.1 or Windows 95, and is 
automated with the exception of user input to a configuration file for setup parameters for each 
program. The input configuration files are updated with a text editor and the various parameters 
are read in during program execution. 

The DPIV acquisition configuration file is called DPIVWRT.CFG and takes the form 
shown below: 


softwareversionnumber 

1.2 

Kodak_exposu relevel ( ms) 

90 

output_datafile_base_name 

datal 

number_of_files 

100 


The configuration file requires three inputs; 1) the Kodak Megaplus exposure level, 2) the output 
datafile base name, and 3) and the number of digital files to capture and store away. The base 
name is used for all the output files with the extension increasing from * 000, * 001, * 002, ...to 
the maximum files to be saved away. 

The DPIV acquisition executable file is called DPIVWRT.EXE and involves 5 basic steps 
(actual code excerpts are shown in italics): 


1) Read in configuration file information: 

readjconfigQ; /* call read_configQ; function call */ 


2) Initialize the Epix frame grabber and Kodak camera: 

if (do _openQ < 0) /* initialize framegrabber and return error if problem V 

return (l); 

pxd_videksetCh'); /* set Kodak camera system to high resolution mode */ 

pxd_videkdo(7,0, exposure); /* capture initializing digital image V 

pxd_video( s', pxdjmbufQ); /* turn frame grabber display off for faster implementation */ 


3) Setup output data file name and begin main loop: 

printf("\nDigitizing/writing fie number:\nW); /* show what file number we are currently saving */ 
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for (q = 0; q < ftlemax ; q++) /* loop through the number of files */ 

{ 

4) Capture 1320x1035 pixel digital image: 

pxdvidekdo(7, 0, exposure); /* capture the digital image to framegrabber */ 

pxdvuleofs’, pxdimbufO ); /* turn display off for faster implementation */ 

onetenhundO; /* set up the string for the output filename */ 

strcpy (filename, ""); 

str cat (filename, base); 

strcaUfilename , "); 

strcat(filenamefilenum2); 

ofp =fopen(fdename, "ab" ); /* get the current filename ready to write */ 

5) Read digital image from framegrabber memory and write to PC memory and then to output file: 

for (j = 0; j < 1035; y++ ) /* loop through the pixel rows one at a time */ 

{ 

yposb = j; /* keep track of current row position */ 

ypose =j + 1; 

pxd_iopen(0, pxd imbufO , 0, yposb , 1320, ypose, Y); /* setup pixel row */ 
pxd_io(0, list, 1320); /* read pixel row V 

fvrite( list, sizeof( unsigned char ), 1320, ofp ); /* write pixel row to file */ 

} 

fclose(ofp); /* close file when done */ 

} 

The data is written in binary form to file in row format at about 7 seconds per 1320x1035 8-bit 
pixel file. 

The DPIV autocorrelation analysis configuration file is called DPIVCOR.CFG and takes 
the form shown below: 


software version number 

1.0 

image_threshold_level(0.9-l. 1) 

1.0 

FP_DC_peak_extent( 1 -20) 

6 

search_box_xmin(0-128) 

45 

search_box_xmax(0- 128) 

110 

search_box_ymin(0-128) 

25 

search_box_ymax(0- 1 28) 

45 

Default_peak_xpos(0-64) 

2 

Dcfault_peak_ypos(0-64) 

2 

Centroid_threshold_level(0-255) 

250 

pixel_scale_factor(microns/pixel) 

1.0 

acquisitionmagnification 

1.0 

laser_pulse_separation(microsec) 

1.0 

input_DPIV_binary_filename 

datal 

outputvectordatafilename 

analyzl 

number_of_files_to_analyze 

10 
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The input configuration file supplied various analysis system parameters needed for 
processing. These included 1) the acquisition magnification factor, 2) image and centroid 
threshold levels, 3) the DC peak correlation extent, 4) correlation search box limits, 5) default 
peak correlation pixels, 6) laser pulse separation timing, and 7) the input and output file base 
names and number of files to analyze. 

The DPIV autocorrelation analysis executable file is called DPIVCOR.EXE and involves 
7 basic steps (actual code excerpts are shown in italics): 

1) Read in configuration file information: 

readjConfigO; /* call readconfigQ; function call */ 

2) Initialize Alacron i860 array processor: 

if (a lope n (DEV) != 0) /* initialize i860 and check for error */ 

errexit ("can't open device %d\n ", DEV); 
aldev (0); 

if (almapload (I860PROGRAM, STACKS1ZE) != 0) /* load analysis program into i860 memory */ 

errexit ("can't load %s\n ”, I860PROGRAM); 


3) Setup output data file name and begin main loop: 
printf("\n \n A nalyzing file number: \n \n "); 

for (qcount = 0; qcount < filemax ; qcount++) /* loop through the number of files */ 

{ 

onetenhundO; /* set up the strings for the output filename */ 

strcpy (filename, ""); /* and input filename */ 

strcatfilename , base); 

str cat (filename , " "); 

strcat(filenamefilenum2); 

strcpy(flename2, ""); 

strcat(flename2, base 2); 

strcat(flename2, ". "); 

strcat(filename2filenun\2); 

ifp = fopenfilename, " rb "); 


4) Read data from file and send to PC memory and to i860 memory: 

read jar ray 0; 

alsetba (VtoP (A bufl), (unsigned char *)buJU, 32000); 
alcall (aladdr (" senddata "), II, A_bufl, A_buf2, A_buJ3, 
A_buf4, A_buf5,A_buf6, A_buJ7, A_buJ8, A_buf9, 
A buflO, I); 
alwait 0; 


/* read data from file */ 
/* send data from PC */ 
/* memory to I860 */ 

/* memory */ 

/* wait for completion */ 


5) Perform FFT and autocorrelation calculations on i860: 

alcall (aladdr (" doff new"), 2, A_bufI28, Ajcountp); /* Run the autocorrelation routine +/ 
alwait Q; /* wait for completion */ 

6) Perform centroid calculation: 

main searchQ; /* Run the centroid extraction routine */ 


1 ) Write data to file: 

fl = fopenftest, "a"); /* set up output filename to write */ 

for (i=0;i<I5;i++) /* loop for 15x19 vector arra v */ 

f 
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/* scale the vectors */ 


for C i-=0;j<J9;j++ ) 

{ 

x [k] = (x[k]*tl0)/(lll*tl2); 
xl[kj = (xl[k]*tlO)/(tl l*tI2); 
x2[k] = (x2[k]*tl0)/(tl I *t 12); 
y[k] = (y[kj*tl0)/(tll*tl2 ); 
yl[kj = (yl[k]*tl0)/(tll*tl2); 
y2[k] = (y2[k]*tl0)/(tl 1 *112); 
countO = (Q *1. 0)+ct[0J); /* keep track of position */ 

count 1 = ((i*1.0)+ct[lj); 
ahsx - countO; 
absy = count 1 ; 

fprintf(fl. "%f%f%f%f%d %f%f%d %f%f%d\n", /* write to file V 

countO, count 1 , x[ k],y[k ],b[k],xl [kj, 
ylfkj, b 1 1 k],x 2 [ k],y 2 [k], b 2[k]); 

k++; 

} 

} 

fclose(fl); /* close output file */ 

The data is written in ASCII format with the top three autocorrelation centroids being saved 
according to peak intensity level which is also saved. A weighted mean centroiding scheme is 
used for sub-pixel centroid extraction. The current autocorrelation algorithm running on the i860 
computes 5 vectors per second which saves away the 15x19 vector array for a single file in 
approximately 1 minute. 
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11 











Figure 9. Ensemble mean velocity all stations with axial component contours. 



Figure 10. Ensemble mean velocity all stations with transverse component contours. 



Figure 1 1 . Ensemble axial velocity variance all stations. 



Figure 1 2. Axial cuts of ensemble axial velocity variance along jet centerline and jet edges. 
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Figure 18. Axial cuts of Turbulence Intensity <v'v’> in Figure 17. 
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Figure 19. Reynold’s Stress <u’v’>. 



Figure 21 . LV, LTA and PiV u-component Standard Deviation 
comparison along jet edge in axial direction. 


Figure 22 LTA and PIV v-component Standard Deviation 
comparison along jet edge in axial direction. 
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